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In recently introduced schematic lattice gas models for vibrated dry granular 



J> . media, we study the dynamical response of the system to small perturbations 
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of shaking amplitudes and its relations with the characteristic fluctuations. 
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Strong off equilibrium features appear and a generalized version of the hue- 



systems and the role of Edwards' compactivity are discussed. 
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Few years ago, Edwards formulated the hypothesis that it is possible to extend to pow- 

O ■ 

ders the methods of standard statistical mechanics an important challenge since powders 
are by definition "non thermal" systems ||. Actually, in granular media the role of a "tem- 
perature" , linked to the concept of Edwards' compactivity |]] , is played by the amplitude of 
external vibrations , and it is possible to properly individuate the corresponding "equilib- 
rium" states However, these systems are typically in off equilibrium configurations, 
as the presence of aging phenomena shows and the above extension of statistical me- 
chanics thus runs across its off equilibrium version. 

A very important issue for both theoretical and practical reasons, is the understanding 
of relations between the response of a system to external perturbations and its characteristic 
fluctuations. Here, we study such a problem in the context of schematic lattice gas models 
|7]-|| recently introduced to describe gently shaken granular systems. We show how in 
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granular media it is possible to formulate a fluctuation dissipation theorem (FDT) where 
the shaking amplitudes play the role of usual temperatures, but in typical off equilibrium 
situations it coincides neither with its usual version at equilibrium nor with the extensions 
valid in off equilibrium thermal systems in the so called "small entropy production" limit 
||10|| . We also discuss the important relations with Edwards' theory. 

The lattice gas model we consider here, the Tetris, schematically describes the effects of 
steric hindrance and geometric frustration on grains in granular materials. Interestingly it 
shows many phenomena typical of granular media under shaking, as logarithmic compaction 
or segregation . It consists of a system of elongated grains, with two possible orienta- 

tions, which occupy the sites of a square lattice tilted by 45°. To avoid overlaps, two particles 
can be nearest neighbor only if they have the right reciprocal orientation. The motion of 
grains in absence of "vibrations" is subject only to gravity so grains only move downwards 
(without overlapping). The effect of vibration is introduced by allowing the particles to 
diffuse upwards with a probability p up and downwards with Pdown = 1 — Pup- The effects of 
geometric frustration on the microscopic dynamics are introduced with a kinetic constraint: 
particles can flip their orientation only if at least three of their neighbors are empty. 

The parameter governing the dynamics is the adimensional quantity x = p up / Pdown, or 
r = 1/ ln(x ) . T plays the same role as the amplitude of the vibrations in real granular 
matter || , and, for not too high amplitudes, a good agreement is found between the model 
and experiments by posing T b ~ a/g with b about 1 (here g is the gravity and a the peak 
acceleration of the shakes in the experiments) f|[3|. In analogy with experiments, in our 
Monte Carlo simulations the shaking of the system is as follows: we prepare the system, 
confined in a box, in a given initial configuration at t = (see below), then we start to 
"shake" it continuously and indefinitely with a given "amplitude" x$. We expect very similar 
results by considering, instead of a single long tap, a series of short taps, as experimentally 
more convenient (see Ref. PJTT||). 
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FIGURES 




FIG. 1. Top left The average grains height ho(t) (filled circles) and h\{t,t w ) (empty squares) 
as a function of time t — t w , of a lattice granular system initially prepared in a uniform fluidized 
state and then shaken at xo = 0.8, and of a replica perturbed after a time t w = 370 by shaking at 
x\ = xq + Axo = 0.802. Top right The "displacement", B(t, t w ) = Co(t, t) — 2Co(i, t w ) + Co(t w ,t w ) 
(where Co(t,t') = (ho(t)ho(t'))), of the above system at xq = 0.8 as a function of t — t w . Bottom 
left The average height ho(t) (filled circles) and hi(t,t w ) (empty squares) as a function of t — t w , 
of a system initially prepared in a static compact state and then shaken at xo = 0.001, and of a 
replica shaken at x\ = .003 after a t w = 333. Bottom right The "displacement", B(t,t w ), of the 
system at xo = 0.001 as a function of t — t w . 

The Tetris can be described with the following Ising Hamiltonian in the limit J — > oo 
(see 0||): (3H = JJ2(ij) fij(Si, Sj)ni7ij — ^ J2t Vi ■ Here g is the gravity constant, yi is the 
z-th particle height, n« = 0, 1 are occupancy variables, Si = ±1 are Ising spin variables for 
the two possible orientations of grains, fij is a function describing the hard core repulsion 
(J = oo): fij(Si,Sj) = if the orientation of neighbors, (Si,Sj), is allowed and = 1 if 
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it is not allowed ( fij(Si, Sj) = 1/2 [SiSj — e^Si + Sj) + 1], where = +1 for bonds along 
one direction of the lattice and e^- = —1 for bonds along the other). The temperature, T, 
of the above Hamiltonian system is related to the ratio x = p up /Pdown via e~ 29//T = x (i.e., 
r = T/g). 

This Hamiltonian mapping shows that in our model the field coupled to grains height 
is the adimensional gravity: g/T = Iti^Xq 1 ^ 2 ). Thus a perturbation to the system, coupled 
to a easy observable, may be introduced by varying xq. We record, with Monte Carlo 
simulations, the dynamical correlation functions and the response of the system to such a 
small perturbation in the "shaking amplitude" . 

Inside a box of fixed size 30 x 60, with periodic boundary conditions along the x-axis 
and rigid walls at bottom and top, we prepare the system in a uniform density initial 
configuration (p = 0.5), corresponding to an highly fluidized state, and then we "shake" it 
at a given amplitude xo- During this process we record the average height h Q (t) = (J2iVi{t)) 
of the grains and the two times correlation function Co(t,t') = (ho(t)ho(t')) , or, the "mean 
square displacement" B(t,t') = Co(t,t) — 2Co(t,t') + Co(t',t') which is a quantity relevant 
to test the FDT theorem. Time is measured in such a way that unity corresponds to an 
average update of all the degrees of freedom in the system, and the statistical averages 
run over 2048 noise and initial configuration realizations. In order to measure the response 
function, we also record the average height of grains, hx(t,t w ), in an identical copy of the 
system (a "replica"), which, after a fixed time t w (typically below we fix t w = 370,3700), 
is shaken at x\ = xq + Axo, i.e., is perturbed by a small increase of the shaking amplitude 
Ax (in what follows Ax = 0.002 @). 

The quantity Ah(t,t w ) = hi(t,t w ) — h (t), i.e., the difference in heights between the 
perturbed and unperturbed systems, is by definition the integrated response. FDT or its 
generalizations concern the relation between Ah(t, t w ) and the displacement, B(t, t w ), which 
is linked to the correlations in the unperturbed system. The simplest version of a possible 
generalization to off equilibrium granular matter of the FDT may be argued from thermal 
systems [fL0| , |l3ll . The integrated response, Ah, should be approximately proportional to B: 
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Ah(t,t w )~ y Air-^Bfat*,) (1) 

Here, A(T _1 ) = T^ 1 — T^f 1 , is the variation between the inverse shaking amplitudes of the 
reference and the perturbed systems and the prefactor, X, is a quantity to be determined, 
in principle function of t w and t themselves. In thermal system, in the limit t,t w —>■ oo, X 
is a piecewise constant, depending just on B(t,t w ) and not on both t and t w flT0| , |i"3 . 



In the specific case of our model we have A(T 1 ) = A(l?) = In ( Xo+Ax ° ) 2 ( and if 



A) 



X = 1 we recover the usual well known equilibrium version of FDT. In the study of glassy 
thermal systems the quantity T/gX =T/X has the meaning of an "effective temperature" 



of the sample, which only for X = 1 coincide with the equilibrium bath temperature |I0 



In the present model, eq. ([!]) seems to be approximately valid, also if in typical off- 
equilibrium situations X slowly depends on t and t w . To outline this, as first we explore the 
high shaking amplitude regime. In Fig.|l| (top left), as a function of t — t w are shown, for 
t w = 370 (analogous results are found for t w = 3700), the height ho(t) of a system "shaken" 
at Xo = 0.8, and h\(t,t w ) of the replica at X\ = Xq + Ax = 0.802. ho and hi decrease in 
time signaling that the system is slowly compactifying while approaching the equilibrium 
density profile. In Fig.|l] (top right), is also shown the "displacement", B(t,t w ), of the same 
unperturbed system. We show in Fig.|2] (top left) the difference Ah(t,t w ) = hi(t,t w ) — h (t) 
as a function of t — t w . In agreement with the expected simple behavior of a gas subject 
to gravity (which reacts to an increase of temperature, or decrease of gravity, by increasing 
its average height), Ah is a positive quantity: the "colder" system (with x = 0.8) actually 
must have a lower equilibrium height respect the "hotter" one (with x\ = 0.8 + Axq). In 
order to check the generalization of the FDT proposed in eq. (HD, in Fig.^| (top right), we 
plot Ah as a function of the displacement B. Actually, Ah is approximately piecewise linear 
in B: after an early transient with X > 1, in the long time regime X is independent on t 
and t w and equal to 1, showing that, in the high "temperature" and low density region, the 
usual equilibrium version of FDT is obeyed. 

In the low Xq region, we already know the presence of strong off equilibrium phenomena 
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@. Thus it is reasonable to expect that the above simple picture drastically changes. In 
FigfJ (middle left), we plot the integrated response, Ah(t,t function of t — t w for 

a system "shaken" at Xq = 0.5 and a replica shaken at Xq + Axq = .502 after a t w — 370 
(analogous results are found for t w = 3700). It is apparent that after an early transient, up 
to the time scales we explored we find a negative response, Ah, in complete contrast with 
the above equilibrium scenario. Actually, both replicas start from a strongly off equilibrium 
state, but while the "colder" is "frozen" and remains longly trapped in metastable states, 
the "hotter" is able to more rapidly escape to approach its asymptote. 

Notice that Ah is negative even after five decades in time and no change in its trend is 
observed. The approach to equilibrium is actually extremely slow ||||: if the region with 
negative response would extend up to infinity (which is hard to say with computer simula- 
tions), this should results in a dynamical breaking of ergodicity introduced by gravity in our 
system of particles interacting just with excluded volume effects. Interestingly these results 



are partially supported by some recent experiments (TTJ, which showed that the asymptotic 
density of a granular system compactified at a low shaking amplitude from a random initial 
configuration, is lower than the analog quantity in a system shaken at a slightly higher am- 
plitude. This result is in contrast with the "equilibrium" measures (the "reversible branch" ) 
in the experiments of Novak et al. ||. The above phenomenon has strong repercussions 
on the FDT. As plotted in Fig|| (middle right), Ah may be approximately plotted as a 
piecewise linear function in B, as from eq. (|l|), but after an initial transient with X > 1, the 
system enters a region with negative X (X ~ —5), corresponding to a negative "effective 
temperature" T/X. The specific value of X, in this region, also slowly depends on t w (see 
below), and we are far from the small entropy production limit known in thermal systems 
where X = X(B). 

A similarly anomalous picture is found in a system shaken at Xq = 0.001 and a replica 
at xo + Axo = .003. To show that the above general results do not depend on the details 
of the initial state, we discus the latter Monte Carlo shaking experiment in a differently 
prepared system, closer to experiments. Now the starting particle configuration is prepared 
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by randomly inserting particles into the box from its top and then letting them fall down, 
with the above dynamics, until the box is half filled. Thus, the system starts from a more 
compact static state. In Fig.|I], the heights of the original system, h , of the replica h\ 
(bottom left), along with the displacement B(t,t w ) (bottom right) are plotted as a function 
of t — t w for t w = 333. This figure shows again a negative response Ah which is plotted as 
a function of t — t w in Fig.|] (bottom left), over six order of magnitudes, for three different 
values of t w (t w = 33, 333, 3330). Three different regions are observed for Ah, but the short 
times transient inflection zone seems to disappear for t w long enough. In agreement with the 
generalized FDT of eq. ([!]), Ah is at long times (10 4 < t — t w < 10 6 ) again approximately 
linear in B, as shown for t w = 333 in Fig.|^ (bottom right) ]14]]. The negative proportionality 
coefficient, X, slowly depends on t w (from X ~ —7 for t w = 33 to X ~ —15 for t w = 3330). 
For systems prepared in the previous fluidized initial state, discussed above, analogous results 
are found. It is important to stress that the above behaviors with negative responses are 
not found if the system starts from an equilibrium initial configuration. 

Interestingly, in the high density or low xq region, we observe an analogous scenario also 
in a similar lattice model, the IFLG |8|], where, in order to describe the motion of grains 
in a disordered environment, quenched disorder is introduced in the Hamiltonian described 
above, but no kinetic constraints are present in the dynamics. 
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FIG. 2. In the figures on the left column, we plot, as a function of t — t w , the average height 
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difference, Ah(t,t w ) = h\(t,t w ) — ho(t), of a reference system shaken at a given xo and a replica 
perturbed after t w by shaking at xq + Axo (Axo = .002). In the right column, in order to check 
the generalized fluctuation dissipation theorem (FDT) of eq. ([j]), we plot the quantity, Ah(t,t w ), 
i.e., the integrated response, as a function of the displacement of the reference system, B(t,t w ). 
The top and middle cases correspond to systems initially prepared in a uniform fluidized state, 
then shaken at different "amplitudes" xq (xq = 0.8 top, xq = 0.5 middle), with replicas perturbed 
after a t w = 370. In the bottom figure, the systems, initially prepared in a static compact state, 
is shaken at xq = 0.001 and its replica is perturbed by Axo = 0.002 after a t w = 33, 333, 3330. As 
an "hot" gas at equilibrium has a higher average height respect to a "colder" gas, Ah should be 
always positive. However, at low xq (xq = 0.5,0.001), negative responses appear. In agreement 
with eq. (|l|), Ah is asymptotically still approximately piecewise linear in B, but, up to the time 
scales we explored, only at xo = 0.8 the equilibrium version of FDT holds. The superimposed 
linear curves indicate the different regimes described in the text. 



As stated, the presence of negative X in the generalized FDT, is a feature of the very 
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far from equilibrium dynamics. Actually, this is not expected for instance in off equilibrium 
thermal systems as glasses or spin glasses, at least in the small entropy production limit 
||10||. Moreover, in spin glasses the parameter X, which is a seemingly dynamical quantity, 



may be asymptotically related to equilibrium static properties from replica theory [10,15 



In granular media, as long as in fragile glasses |16|], the static correspondent of X, if any, 
is still missing. It is useful to underline that while in glasses an experimental measure of 
the generalized FDT may be non trivial, in granular media, as shown above, this should be 
reasonably simple. 

The interesting finding in our models for granular matter of a negative X, corresponding 
to negative responses and "effective temperatures" , may be clarified by stressing the relations 
with Edwards' theory |L]. 

Actually, the present approach to granular media based on a standard Hamiltonian 
formalism, is very close in spirit to the more general statistical mechanics approach proposed 
by Edwards |1| . The hard core repulsion term in the above Hamiltonian H has a role similar 
to Edwards' volume function W (where a constraint of mechanical stability is explicitly 
present). The fundamental control parameter in Edwards' theory which corresponds to our 
adimensional temperature, T/g, is the compactivity AX^. Edwards and Grinev || guess that 
the compactivity is related to the experimental "temperature" of a shaken granular medium 
according to the following "fluctuation-dissipation" relation: XX e ~ {a/g) 2 (as above, a is 
the shake and g is the gravity acceleration). This statement is the analog of our suggestion, 
{T/g) b ~ a/g, which is based on the comparison of Monte Carlo and experimental data. 

In Ref. [|J is proposed that the state of a shaken granular medium in the equilibrium 
regime, corresponding for instance to configurations on the so called "reversible branch" in 
the experiments of Ref. ||, is characterized by a positive compactivity, and, by extrapolation, 
the off equilibrium dynamics, corresponding to the experimental "irreversible branch" of 
Ref. H, has a negative compactivity, Xe < 0. This observation is in agreement with our 
discover of different regions with positive as long as negative effective temperatures, T/X, in 
the study of the generalized FDT in granular media. A negative off equilibrium compactivity 
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may correspond to the necessity of the system to cross states with "higher entropy" in 
order to lower its volume under shaking. Actually, the measure of integrated responses v.s. 
correlation functions in granular materials may open the way to a direct experimental access 
to Edwards' compactivity for a clear settlement of a statistical mechanics for such systems. 

In summary, in order to understand the off equilibrium statistical mechanics of powders, 
we have studied, in schematic lattice gas models for vibrated dry granular media, the dy- 
namical response functions to small shaking amplitude perturbations and their relations to 
characteristic dynamical fluctuations. Strong off equilibrium features appeared, as long time 
regions with negative response functions, different from those observed in the small entropy 
production limit of thermal glassy systems |10| , along with the necessity to introduce a gen- 
eralized version of the fluctuation dissipation theorem. The novel properties we have found, 
as the presence of negative effective temperatures, are consistent with Edwards' theory of 
powders and demand important experimental check. 

We thank INFM-CINECA for computer time on Cray-T3D/E. Work partially supported 
by the TMR Network Contract ERBFMRXCT980183 and MURST-PRIN 97. 
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